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TURBULENCE IN GALAXY CLUSTER CORES: A KEY TO CLUSTER BIMODALITY? 



Astronomy Department 



o 

(N 



(N 

43 
6 

CO 

> 

in 
oo 

(N 

On 
O 



13 



Ian J. Parrish 1 , Eliot Quataert, and Prateek Sharma 1 

Theoretical Astrophysics Center, 601 Campbell Hall, University of California Berkeley, CA 94720, USA; 
iparrish@astro.berkeley.edu 
Draft version March 30, 2010 

ABSTRACT 

We study the effects of externally imposed turbulence on the thermal properties of galaxy cluster 
cores, using three-dimensional numerical simulations including magnetic fields, anisotropic thermal 
conduction, and radiative cooling. The imposed "stirring" crudely approximates the effects of galactic 
wakes, waves generated by galaxies moving through the intracluster medium (ICM), and/or turbulence 
produced by a central active galactic nucleus. The simulated clusters exhibit a strong bimodality. 
Modest levels of turbulence, ~ 100 km s _1 ~ 10% of the sound speed, suppress the heat-flux-driven 
buoyancy instability (HBI), resulting in an isotropically tangled magnetic field and a quasi-stable, 
high entropy, thermal equilibrium with no cooling catastrophe. Thermal conduction dominates the 
heating of the cluster core, but turbulent mixing is critical because it suppresses the HBI and (to a 
lesser extent) the thermal instability. Lower levels of turbulent mixing (< 100 km s _1 ) are insufficient 
to suppress the HBI, rapidly leading to a thermal runaway and a cool-core cluster. Remarkably, then, 
small fluctuations in the level of turbulence in galaxy cluster cores can initiate transitions between 
cool-core (low entropy) and non cool-core (high entropy) states. 

Subject headings: convection — galaxies: clusters: general — instabilities — MHD — plasmas 



1. INTRODUCTION 

The cooling time in the intracluster medium (ICM) 
of galaxy clusters is often < 0.1 -1 Gyr near the cen- 
ter of the cluster (Sarazin 1986). X-ray spectroscopy 
shows, however, that the majority of the plasma is not 
in fact cooling to tem peratures well below the m ean clus- 
ter temperature (e.g.. iPeterson fc Fabianl 120061) . Under- 
standing the processes responsible for heating and sta- 
bilizing cluster plasmas is a central problem in galaxy 
formation. Some of the most promising energy sources 
include a central acti ve galactic nucleus (AGN) (e.g., 
iBinney fc Taborl I1995D. thermal conduc tion from large 
radii (e.g.. iNaravan fc Medvedevl |2001[ ). or dynamical 
friction and/or turbulence generated by the motion o f 
substructure through the cluster (e.g.. iKim "elall[2005h . 
In this Letter we show that although neither of the latter 
two mechanisms works on its own, together they have 
novel implications for the thermodynamics of the ICM. 

The central parts of clusters are unstable to a con- 
vective instability driven by the anisotropi c flow of 
heat along magnetic field lines (the HBI; lOuataertl 
2008). Simulations of the HBI show that it satu- 
rates by preferentially reorienting the magnetic field 
lines to be perpendicular to the temperature gradi- 
ent, thus r educing the effective radial conductivity of 
the p l asma dParrish fc Quat aert 2008; Bogda novic et all 
I2009t iParrish et alj I2009L hereafter PQS). This exacer- 
bates the cooling catastrophe by making it difficult to 
tap into the thermal bath at large radii, particularly for 
clusters with low central entropies and short central cool- 
ing times. 

Previous simulations of the HBI in clusters have fo- 
cused on idealized problems i n which the HBI was the 
primary source of turbulence (|Sharma et al.1 l2009al also 
studied the role of convection driven by cosmic-rays). 
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There are, however, many additional sources of turbu- 
lence in clusters, including major mergers, the motion of 
galaxies through the ICM (galaxy wakes), and AGN jets 
and bubbles. Cosmological hydrodynamic simulations 
of cluster formation find that turbulence can contribute 
~ 1-10% of the total pressure even in relaxed clusters, 
with the turbulent pre ssure declining a t small radii to- 
wards the cluster core (|Lau et al.ll2009f h 

In this Letter, we consider a simple model for the in- 
terplay between turbulence, anisotropic thermal conduc- 
tion, and radiative cooling in galaxy cluster cores: we 
externally "stir" the ICM in our previous global cluster 
core simulations (e.g., PQS) in order to mimic the ef- 
fects of the various sources of turbulence noted above. 
The limitations of this approach are discussed in jjH 

Nea r the completion of this work, iRuszkowski fc Ohl 
(2009) presented results similar to those found here using 
independent numerical techniques and cluster models. 

2. METHODS 

We solve the equations of MHD using the Athen a 
MHD code (|Gardiner fc Stond 12001 [Stone et alJl200l . 
with the addition of anisotropic thermal condu ction 
(jParrish fc Stond [20051 : ISh anna fc Hammettl I2007T ) and 
optically thin cooling (see eqs [8]-[12] of PQS). In par- 
ticular, the conductive heat flux is given by Q — —Ks p bb- 
VT, where ks p is the Spitzer thermal conductivity and 

b is a unit vector along the magnetic field. We use the 
iTozzi fc Normanl (|200lh cooling curve and a temperature 
floor of T = 0.05 keV, below which UV lines become im- 
portant. 

Our initial condition is a clus ter inspired to resembl e 
that of Abell 2199 as observed in lJohnstone et al. (2002). 
We use a static NFW gravitational potential with a scale 
radius of r s = 390 kpc and a mass of Mq = 3.8 x 10 14 M Q . 
The simulations are carried out on a Cartesian grid in a 
computational domain that extends from the center of 
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the cluster out to 240 kpc. The simulations are (128 3 ), 
corresponding to a resolution of 3.4 kpc. Resolution 
studies indicate that our conclusions are not sensitive 
to resolution. We begin with an ICM that is in both 
hydrostatic and thermal equilibrium, with conduction 
balancing cooling; simulations without initial thermal 
equilibrium showed large thermal transients. The result- 
ing model cluster has a central temperature and electron 
density of ~ 1.3 keV and ~ 0.04 cm -3 , respectively, and 
a temperature and density of 5 keV and 2 x 10 -3 cm~ 3 
near 200 kpc. The magnetic field is initially tangled, 
with (\B\) = 10~ 9 G and a Kolmogorov power spectrum. 
Further details of our initial conditions can be found in 
sections 3-4 of PQS. 

Our model clusters do not include the turbulence 
that would realistically be generated by the hierarchi- 
cal growth of structure or a central AGN. To crudely 
account for this, we continuously add an additional ran- 
dom velocity to our cluster models. We drive the veloc- 
ity r ^elds_Jn_^burierspa£e_^ing the methods described 
in iLemaster fc Stonel ([2009) ■ For a given driving length 
scale L, we drive velocities with a flat spectrum in Fourier 
space on scales corresponding to L ± 10 kpc. We clean 
the spectrum so that the motions are incompressible and 
Fourier transform to real space normalizing to the desired 
level of turbulence. This results in rms velocities, Sv, that 
are uniform throughout the cluster. Note that our stir- 
ring is statistically steady in both time and space. This 
is not necessarily a good approximation in clusters — we 
return to this point in []4l 

Our fiducial turbulence parameters are Sv ~ 50- 
100 km s _1 and L ~ 40-100 kpc. This corresponds to 
turbulence that contributes a few % of the total pres- 
sure in the cluster core. The strength of the turbulence 
produced by galaxies moving through the ICM can be 
estimated by calculating the total po wer dissipated b y 
dynamical friction (e.g., eqn. [4] of lKim et all l2005h . 
Assuming that this energy is ultimately dissipated via a 
turbulent cascade, we find 
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where M g is the mass of a typical galaxy, N g is the num- 
ber of galaxies, R is the size of the region of interest, 
and c s is the sound speed of the ICM. For five 10 n Af Q 
galaxies within 200 kpc and a turbulent scale of L ~ 40 
kpc, equation ([1]) implies Sv ~ 0.08c s , consistent with 
our fiducial numbers quoted above. In reality, the gener- 
ation of turbulence by galaxies moving through the ICM 
will be more subtle, with some of the energy going into 
sound waves and gravity waves, and some being confined 
to galactic wakes correlated with the orbits of galaxies. 

3. RESULTS 

To illustrate the effect that turbulence can have on the 
thermal evolution of galaxy clusters, Figure [T] shows the 
late-time azimuthally averaged temperature profile for 
simulations of the same cluster with different rms tur- 
bulent velocities, Sv, and thus also different turbulent 
heating rates eturb — p(Sv) 3 /L. For low turbulent veloc- 
ities, the evolution is similar to that found previously by 
PQS: the HBI reorients the magnetic field to be perpen- 
dicular to the radial temperature gradient, shutting off 
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Fig. 1. — Azimuthally-averaged temperature profiles for iden- 
tical cluster cores with different imposed rms turbulent veloci- 
ties (see legend); the driving scale is fixed at L = 40 kpc. A 
strong bimodality in cluster properties results. Stronger driving 
(Sv > fOO km s" 1 ) leads to a roughly stable thermal profile (shown 
at fO Gyr). Somewhat weaker driving leads to a cooling catas- 
trophe (shown at the onset of the cooling catastrophe). Eqn. |21 
quantifies the level of turbulent mixing required to suppress the 
HBI and transition from the low temperature, low entropy state, 
to the high temperature, high entropy state. 



heat conduction from large radii and thus precipitating a 
cooling catastrophe in the cluster core. The cluster core 
reaches the temperature floor at ~ 1.9 Gyr in our lowest 
Sv simulation. For larger Sv > 100 km s _1 ~ 0.1 c s , how- 
ever, the dynamics changes completely. The magnetic 
field remains relatively isotropic at all times, indicating 
that the HBI is no longer acting effectively. Moreover, 
the cluster reaches a stable equilibrium, with the tem- 
perature profiles shown in Figure [T] remaining roughly 
the same for the last ~ 5-7 Gyr of the simulation. It is 
important to stress that even for Sv ~ 100 km s _1 , the 
heating rate due to the turbulence is negligible compared 
to the cooling rate and thus the turbulence is energeti- 
cally unimportant for the thermal properties of the clus- 
ter (0. Note also that the central temperature increases 
slightly as the turbulent energy increases in Figure [T] we 
attribute this to the increased advective (turbulent) heat 
transport associated with the higher turbulent velocities. 

To quantify in more detail the effect of turbulence on 
the evolution of cluster plasmas, Figure [2] shows the tem- 
perature and magnetic field direction as a function of ra- 
dius at several different times for two simulations (labeled 
A and B) whose properties are summarized in Table Q] 
(case C in Table Q] is discussed below) . The simulations 
are again of identical clusters and both include turbu- 
lence with Sv ~ 100 km s _1 , near the threshold for the 
transition from stability to instability in Figure [TJ The 
simulations differ in that case A has a turbulent correla- 
tion length of L = 40 kpc, while B has L = 100 kpc; as 
a result, simulation A has a heating rate eturb that is 2.5 
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TABLE 1 

Properties of the Fiducial Simulations' 1 



Run 


(eturb) (erg cm 3 s 1 ) 


L (kpc) 


Sv (km s 1 ) 


tcool (Myr) 


*HBI (Myr) 


toddy (L) (Myr) 


Wdy (Af) (Myr) 


K§ nal (keVcm 2 ) 


A 


7.5 x icr 30 


40 


115 


400 


100 


400 


360 


53 


B 


3.0 x icr 30 


100 


112 


400 


100 


1000 


490 


0.06 


C 


3.0 x 10" 29 


100 


250 


400 


100 


390 


195 


110 



a Timescales are estimated at the initial time, while the turbulent velocity and central entropy are measured in the saturated state for runs A & 
C, and just before the cooling catastrophe for run B. The Field length is estimated near the cluster center, while e tur b is volume averaged. 
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Fig. 2. — Azimuthally-averaged profiles for simulations with 
imposed turbulence having the same Sv ~ 112 km s — 1 , but dif- 
ferent correlation lengths. Run A (left) has L = 40 kpc, while 
run B (right) has L = 100 kpc (see Table [TJ Top: Temperature 
profiles. Bottom: Magnetic field angle relative to the radial: 60° 
corresponds an isotropically tangled magnetic field and an effective 
radial conductivity ~ 1/3 Spitzer. Run A (left), with the shorter 
turbulent mixing time, reaches a stable state averting the cool- 
ing catastrophe. In Run B (right), the HBI persists, shutting off 
conduction and initiating a cooling catastrophe. 

times higher and an eddy turnover time on scale L that 
is 2.5 times shorter. 

Despite their similarities, the top panels of Figure [2] 
show dramatic differences in the evolution of the clus- 
ters' radial temperature profile. Simulation B proceeds 
with little regard for the turbulent driving. The HBI 
reorients the magnetic field, reducing the effective ra- 
dial thermal conductivity, and hastening the onset of the 
cooling catastrophe at ~ 2.2 Gyr. On the other hand, 
simulation A is dramatically affected by the turbulent 
driving, just like the high Sv simulations in Figure [TJ 
In this case the turbulence effectively shuts off the HBI, 
with the mean angle of the magnetic field from radial 
fluctuating about the isotropic value of (Ob) ~ 60°. 

In addition to suppressing the effects of the HBI, the 
presence of "sufficient" turbulence also appears to mod- 
ify the thermal stability of the cluster. Case A in Fig- 
ure [5] shows that the cluster reaches a statistically stable 
thermal equilibrium that survives for longer than the age 



of the universe. The stability of the new equilibrium is 
illustrated by the fact that the cluster's central temper- 
ature undergoes slight oscillations about a new equilib- 
rium value. 

4. INTERPRETATION AND DISCUSSION 

It is important to stress that for all of the simulations 
presented here, the heating produced by the turbulence 
is energetically subdominant (this is also consistent with 
the modest dynamical frictio n heating in cluste rs inferred 
using observed galaxies; e.g.. lKim et al.l l2005). For sim- 
ulation A in Figure [21 e.g., the turbulent heating rate per 
unit volume at the center of the computational domain 
is < 2% of the central cooling rate throughout the evolu- 
tion of the cluster. As a result, the physics important for 
the results in Figures [TJ & [2] includes anisotropic thermal 
conducti on, the HBI , the thermal instability of cluster 
plasmas ([Field! Il965l ) . and the mixing produced by the 
turbulence — but not the heating by such turbulence. 

In general the global thermal instability of cluster plas- 
mas in the presence of thermal conduction manifests 
itself as either catastrophic cooling in the core of the 
cluster, or overheating and th e approach to an isother- 
mal temperature profile (e.g., IConrov fc Ostrikerl 120081 : 
PQS). Which of these is realized in a given problem 
depends in part on the initial state of the system and 
boundary conditions. In cluster simulations without ex- 
ternally imposed turbulence, the HBI biases the nonlin- 
ear evolution of the thermal instability towards the cool- 
ing catastrophe by thermally decoupling the core from 
larger radii. 

Our simulations show that if turbulence is sufficiently 
strong in cluster cores it can effectively shut off the 
HBI, leaving the magnetic field tangled and relatively 
isotropic, and the cluster with a quasi-steady, not-quite- 
isothermal temperature profile (Figs. [TJ&[5]). Quantita- 
tively, turbulence with Sv > 100 km s _1 or a Mach num- 
ber > 0.1 appears sufficient. More precisely, we believe 
that the critical criterion is ([Sharma et al.ll2009bl ) 

teddy (L) ~ — < £<HBI ~ £ ( g—j— J (2) 

where £ e ddy(-^) is the timescale for the turbulence to mix 
the plasma at the outer scale, £hbi is the HBI growth 
time, g is the local gravitational acceleration in the clus- 
ter, and £ is a dimensionless constant that must be de- 
termined from the simulations. Note that because the 
mixing timescale is smaller on smaller scales in a Kol- 
mogorov cascade, the timescale inequality in equation 
$2$ is the most difficult to satisfy at the outer scale. 

Figure [5] demonstrates explicitly that a given value of 
Sv is not in fact sufficient to suppress the HBI and halt 
the cooling catastrophe. Rather, this only occurs in run 
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A, which has a smaller correlation length L and shorter 
mixing time than run B (Table Q}. Table [T] shows the 
properties of a third simulation not in Figure [U in which 
the correlation length is the same as in run B (L = 100 
kpc), but the turbulent velocity is larger. This combina- 
tion again satisfies equation ([2]) and so the evolution is 
qualitatively similar to run A. Because the HBI growth 
time is ~ 100 Myr at ~ 50 kpc in these models, our 
numerical results correspond to £ ~ 5-8 in eq. ([2]). 

In addition to its effects on the HBI, turbulent mix- 
ing can also significantly modify the thermal stability 
of cluster plasmas. Independent of turbulent mixing, 
thermal instability is stabilized on small scales (along 
the magnetic field) by thermal conduction. The critical 
length-scale below which conduction suppresses thermal 
instability (the Field length) is given by 



A, 



4tt 2 Tk(T) 



n e n p A(T) 



1/2 



53 kpc 



K 



15 keVcm 2 



3/2 



(3) 



where A(T) is the cooling function and we have used the 
full Spitzer conductivity in the final expression; in the 
second equality we have also assumed for simplicity that 
the cooling is pure bremsstrahlung so that Xf can be 



Tn, 



-2/3 



expressed solely in t erms of the entropy K 
([Donahue et aiJl2005[ ). 

Absent turbulence, fluctuations with length-scales > 
Xp are unstable on a cooling time t coo \. If, however, 
turbulent mixing is sufficiently rapid, i.e., if 

teddy (Xf) J$ icool, (4) 

then turbulent mixing will suppress the thermal insta- 
bility on the scale of the Field length, although larger 
length-scale fluctuations may remain unstable; local sim- 
ulations of thermal instability in the presence of back- 
ground turbulence confirm this intuition (these will be 
presented elsewhere) . Table Q] shows that for both runs 
A and B in Figure [2j t e ddy(Xp) ~ i C oolj that is, the 
modest turbulence levels considered here are capable of 
significantly changing the dynamics of the thermal in- 
stability in cluster plasmas. One extreme limit of this 
is the possibility that turbulent mixing of hot gas from 
large radii with cooler gas fro m small radii can help pre- 
vent the cooling catastrophe ([ZuHone et al.l [20091) . In 
our stable simulations, however, (e.g., Case A) this is 
not realized: thermal conduction (not turbulent mixing) 
provides the dominant source of heating at small radii. 
The key role of the turbulent mixing is that it suppresses 
the HBI, isotropizes the magnetic field, and helps sup- 
press the thermal instability by mixing the plasma before 
it can cool. This dynamics is qualitatively analogous to 
the critical role that turbulence plays in mixing and dis- 
rupting laminar conductive flames in the combustion and 
Type la supernova contexts (e.g.. iPetersl 2000; W ooslevl 
120071 respectively). 

In our simulations, the interaction between turbulence, 
the HBI, and cooling leads to a strong bimodality in the 
cluster properties (e.g., temperature profiles). Figure [1] 
shows that runs with moderately strong turbulence (sat- 
isfying eqn. [2] ) reach a quasi-stable thermal equilibrium 
averting the cooling catastrophe. By contrast, runs with 
slightly weaker turbulence — Sv smaller by just ~ 25% — 
progress to a cooling catastrophe on a timescale as short 
as a few central cooling times (much like the p ure HBI 
simulations of PQS and iBogdanovic et al.|[2009l) . 



It is tempting to relate this behavior to the observed 
variety in galaxy cluster properties. Observationally, 
clusters show a bimodality in their central gas entropies 
and cooling times, with lower entropy clusters preferen- 
tially having more sta r formation and Ha emission, and 
more powerful AGN (jVoit et al.1 120081 : ICavagnolo et al.1 
2009); the transition occurs at ~ 30 keVcm 2 . This bi- 
modality is closely related to the well-known fact that 
clusters c ome in both cool co re and non cool core vari- 
eties (see lHudson et al.[|2009l) . 

The observed b imodality in cluste r properties is not 
well understood. iBurns et al.1 (|2008h argued that early 
major mergers could pre yent the formation of cool core 
clusters. Alternatively, iGuo et ail (|2008f) showed us- 
ing both ID time dependent models and a global sta- 
bility analysis that the combination of AGN feedback 
and scalar conduction can produce a bimodal population 
of stable cluster models; in their models AGN heating 
largely balances cooling in lower entropy clusters while 
conduction is more important in higher entropy clusters. 
They further suggested that AGN feed back could trans i- 
tion clusters from low to high entropy (|Guo fc Ohll2009h . 

Our simulations demonstrate that a low entropy cool- 
core cluster can transition to a significantly higher en- 
tropy state: run A initially has a central entropy of 
Kq ?s 11 keVcm 2 while its final central entropy is 
Kq w 57 keV cm 2 . Run C is the same as run B but with a 
higher Sv; its final central entropy is Kq ~ 110 keVcm 2 . 
This increase in central entropy is a consequence of run- 
away conductive heating at the roughly fixed pressure 
required for hydrostatic equilibrium. Physically, such a 
transition could occur if a cluster initially had little tur- 
bulence and inefficient conduction (because of the HBI) , 
but was stabilized by a central AGN. The sudden onset 
of turbulence satisfying equation ([5]) — produced by an 
infalling galaxy or the AGN — would isotropize the mag- 
netic field and suppress the HBI. The cluster would then 
evolve as in Figure [5] (left panel; case A) to a higher en- 
tropy state. If the turbulence in the cluster core later 
died away, the HBI would rapidly rearrange the mag- 
netic field (Fig. [TJ, leading to cooling of the core after 
~ 1 Gyr and (by assumption) increased AGN activity 
that would again stabilize the cluster in a cool-core state. 
This demonstrates that modest levels of turbulence in 
cluster cores (Sv ~ 100 km s _1 ; eqn. [5]) can have a 
dramatic affect on their thermal evolution. Note also 
that the energy required to generate turbulence capable 
of suppressing the HBI and thermal instability is far less 
than that req uired to directly increase the entropy of the 
cluster (as in lGuo fc Ohl[2009h . 

Our calculations are based on an overly-simplified 
treatment of turbulence in galaxy cluster plasmas. Real 
turbulence in clusters is likely to be more intermittent 
in space and time than our model (31]), depending on, 
e.g., the distance to an AGN jet or a galactic wake. 
Waves generated by AGN jets and bubbles and/or galax- 
ies moving through the ICM may produce reasonably 
volume-filling turbulence in cluster cores, but this needs 
to be studied in detail. Provided that the turbulence 
is replenished on i coo i ~ 0.1-1 Gyr, our results will be 
relatively unchanged. Otherwise, the HBI and thermal 
instability will proceed unchecked. This temporal and 
spatial intermittency of turbulence in cluster cores may 
ultimately prove to be a positive feature, not a "bug," 
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of our model: as described above, modest changes in the 
level of turbulence in clusters can produce rapid and dra- 
matic changes in the thermal structure and stability of 
the cluster, to the point of initiating transitions from low 
to high entropy states (and vice- versa) . Overall, the sub- 
tle interaction between turbulence, the HBI, and cooling 
in galaxy cluster cores has a surprisingly large impact 
on the thermal properties of the ICM. The critical role 
of the turbulence is not the small amount of turbulent 
energy dissipated (which is <C the cooling luminosity); 
rather, it is the fact that turbulent mixing can suppress 
both the HBI and the thermal instability in cluster cores 
(see eqs. [2] and 0], respectively). 
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